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^ : Abstract 
O 

' The depinning of a domain wall in ferromagentic metal via macroscopic 



quantum tunneling is studied based on the Hubbard model. The dynamics of 
' the magnetization verctor is shown to be governed by an effective action of 

Heisenberg model with a term non-local in time that describes the dissipation 
^ ' due to the conduction electron. Due to the existence of the Fermi surface 

there exists Ohmic dissipation even at zero temperature, which is crucially 
^ ' different from the case of the insulator. Taking into account the effect of 

pinning and the external magnetic field the action is rewritten in terms of a 
collective coordinate, the position of the wall, Q. The tunneling rate for Q is 
calculated by use of the instanton method. It is found that the reduction of 
the tunneling rate due to the dissipation is very large for a thin domain wall 
with thickness of a few times the lattice spacing, but is negligible for a thick 
domain wall. Dissipation due to eddy current is shown to be negligible for a 
wall of mesoscopic size. 
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I. INTRODUCTION 



Domain wall in a ferromagnet is a soliton connecting two stable spin configurations 
separated by the energy barrier due to the anisotropy (Fig. |l|). It is generally pinned by 
defects where the spins are easier to rotate. 
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FIGURES 



FIG. 1. A configuration of magnetization of a planar domain wall (a) where the spin configu- 
ration spatially varies in x-direction, which is perpendicular to the easy plane of the spin. A planer 
wall with the easy axis perpendicular to the wall is shown in (b) [0. 

The wall gets depinned by the external magnetic field (Fig. This depinning process 
involves the motion of an extended object, but it can be simplified if one neglects the bending 
of the wall and describes the process as a motion of a collective varialbe, the position of the 
wall, Q. The potential for Q created by defects ( Fig. |^(a)) is deformed by the external 
magnetic field if, and the position of the wall at the pinning center becomes metastable 
(Fig. §(b)). The barrier height is reduced as the field is increased and finally vanishes at 
the coercive field, H^. At high temperature the depinning occurs by the thermal fluctuation 
of the variable Q overcoming the barrier, and at low temperature it is due to the quantum 
tunneling of Q. Since the wall contains a large number of spins it is a macroscopic quantum 
tunneling (MQT). 

This MQT of domain wall was investigated theoretically by Stamp IQ-Q. His considera- 
tion on an insulating magnet in three dimensions is based on the ferromagnetic Heisenberg 
model with a uniaxial anisotropy. By use of a classical solution of planar domain wall, 
the hamiltonian describing the position of the wall is obtained. The pinning potential by 
a defect is assumed to be of a short range. The experiments of the tunneling of the wall 
from the metastable minima are carried out in a magnetic field close to this coercive field 
{H < He) for the rate to be large enough to be observed. Stamp obtained the expression of 
the tunneling rate as functions of the coercive field, saturation magnetization, the number of 
spins in the wall and the external field, and concluded that depinning of a wall due to MQT 
is observable even for a large wall containing about 10^° spins. For observing the MQT, 
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highly anisotropic materials (such as those containing rare-earth) would be suitable, since 
the wall will be small (thin) in these materials and also the crossover temperature from the 
thermal to the quantum regime can be higher. 



FIG. 2. (a):A domain wall pinnied at a defect. (b):Depinning of the wall by the external 
magnetic field. 



FIG. 3. The potential for the domain wall in terms of Q, the position of the center of the wall; 
(a) without the magnetic field and (b) with the magnetic field. 

In macroscopic quantum systems, the coupling to the environment, which results in 
dissipation, may become significant, and the effect needs to be estimated. In their seminal 
paper Caldeira and Leggett presented a formalism that can incorporate the dissipation 
in quantum mechanics by use of the imaginary time path integral. There the dissipation is 
expressed by the action non-local in time. Based on this formulation, Caldeira and Leggett 
investigated the quantum tunneling of a macroscopic variable, and found that the dissipation 
generally reduces the tunneling rate. 

As sources of dissipation in the MQT of domain wall. Stamp considered the magnon (spin 
wave) and phonon. Because of the gap due to the anisotropy, the Ohmic dissipation due to 
the spin wave vanishes at absolute zero and is quantitatively negligible at low temperature. 
The phonon was also shown to be irrelevant. The effects of conduction electrons were briefly 
touched by Chudnovsky et al. on a phenomenological ground. Concerning the dissipation 
in single domain magnets, the magnetoelastic coupling to phonons was considered, but the 
effect turned out be negligible in actual situations . The effects of spin wave have not been 
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discussed, but are expected to be negligible since this degree of freedom is frozen out at very 
low temperature due to the anisotropy gap. Recently it was claimed that the dissipation 
due to hyperfine coupling to the nuclear spin is significant [0. 

The depinning of the domain wall is observed by measuring the relaxation of the mag- 
netization, and MQT is indicated by the relaxation rate going to a finite constant value as 
T — > 0. The first indication of MQT was obtained in a bulk ferromagnet of SmCo3.5Cui.5 
P|, but detailed study was carried out on small ferromagnetic particles of Tbo.5Ceo.5Fe2 of 
diameter of about 150A f^. Small particles are suitable for such experiments since in the 
bulk sample many walls with different size participate in the relaxation process. The tem- 
perature independent relaxation was observed for T < 0.6K, and the result was claimed to 
be consistent with theory of MQT of a domain wall in a ferromagnet without dissipation. 

Although the materials employed in these experiments are generally matallic, there has 
been no theoretical work paying full attention to this fact. In metals the behavior of the 
magnetic object may not be so simple as in insulators. Actually the magnetization arises 
due to the polarization of the electron and at the same time these electrons have dissipative 
effects for the magnetization vector. To take into account these points in a systematic 
way our analysis |TD[ of the MQT of a domain wall is based on the Hubbard model. We 



consider the case of absolute zero, T = , since we are interested only in the quantum 



tunneling present at low temperature. The calculation is carried out in the continum |]1T 
The magnetization vector is expressed as an expectation value of the electron spin operator. 
For the treatment of a slowly varying field like a domain wall, the locally rotated frame for 
the electron is convenient. By use of this frame, the effective action that describes the low 
energy dynamics of magnetization vector is obtained by integrating out the electron degrees 
of freedom in the imaginary time path integral. As far as low energy behavior is concerned, 
the modulation of the magnitude of the magnetization is irrelevant, and only its direction 
can be a dynamical variable. Then the part of the effective action that is local in time {i.e., 
instantaneous) has the same form as that of a ferromagnetic Heisenberg model; i.e., there is 
no formal difference from the case of the insulator but the parameters have different physical 



origin. On the other hand, the non-local {i.e., retarded) part of the effective action, which 
describes the dissipative effect of itinerant electron on the motion of the magnetization, is 
crucially different from the case of an insulator. Due to the Stoner excitaion, which is the 
gapless excitation of spin flip across the fermi surface, the Ohmic dissipation is present even 
at zero temperature. It is shown that this dissipation reduces the tunneling rate very much 
particularly for a thin domain wall of thickness comparable to the inverse of the difference 
of the fermi momenta {kpi^ — A;fj,)~^- In the case of strong ferromagnet, on the other hand, 
dissipation resulting from the Stoner excitation is reduced. In such a case the existence of 
non- magnetic band in actual systems is to be taken into account. In contrast to these effects 
from the Stoner excitation, the effects of Ohmic dissipation due to charge current (eddy 
current) will be found to be negligible in a meso- or microscopic wall we are interested in. 

This paper is organized as follows. §|I| is devoted to the derivation of the effective 
action for the magnetization vector on the basis of the Hubbard model. In ^JJ\ the MQT 
of a domain wall is studied based on this action. The expression of the tunneling rate 
including the dissipative effect from itinerant electron is obtained there. The dissipative 
effects resulting from the non-magnetic band are calculated in § [rV| , and in §|V| the effects 



due to the eddy current are estimated. Discussions and conclusion are given in § [VT| . 

II. DERIVATION OF AN EFFECTIVE ACTION 

A. Model 

We consider the Hubbard model in the imagnary time path integral whose Lagrangean 
is given by |]T2 



where Cxo- is an electron operator at site x with spin a{= ±) and rixo- = cl_cr'^yia, ep being 
the chemical potential. For simplicity, a parabolic dispersion with mass m is assumed for 
band energy. The ferromagnetism is induced by the Coulomb repulsion (f/) term, which is 
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rewritten by introducing the magnetization vector M by use of the Hubbard- Stratonovich 
transformation. The vector is expressed as M(x, r) = M(x, r)n(x, r), where n is a slowly 
varying unit vector that represents the direction of the magnetiazation and the magnitude 
of the magnetization is given by 

M(x, r) =< (cVc) >x ■n(x, r). (2) 

The Lagrangian is then written as 

L = Y: clidr + ek)ck. - f/ E M(x)(ct^c). + ^ ^ Mi^)'. (3) 

The spatial variation of n(x) accompanying with a domain wall is much slower compared 
to the inverse fermi momentum of the electron kp^. For the analysis of such a slowly varying 



field, the local frame of electron |]13[ is convenient for the perturbative treatment of the 
interaction between the domain wall and the conduction electron to be valid. In this frame 
the 2;-axis of the electron is chosen in the direction of the local magnetization vector n(x). 
The electron operator in the new frame axo- is related to c^a in the origial one as 
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a^a- = o- COS -Cx,a + e~""^ sin ^Cy,-a (4) 

where {6,(j)) are the polar coordinates of vector n(x, r). We neglect in the following calcu- 
lation the spatial variation of the magnitude M(x) of the magnetization [i.e., M(x) = M), 
since the fluctuation of M(x) has a finite mass. Thus the relevant degree of freedom at low 
energy is only the variation of the angle, {6, 0), of the magnetization. The a^a electron is then 
polarised uniformly in the 2;-direction, and it interacts with space-time variation of the mag- 
netization vector. The interaction term arises from the kinetic energy term c^c+ 1 Vc^/ (2m), 
and is written as 

Hint = E f-<^(l - COS e){a^a,a) + ^ - (±^ - ^ sin ^0)e^^'^(aV±a) 
X 1-2 _i_ 2 

+ 4^Q(V^)' + (l-cos^)(V0f) (ata) 



+ -cos^)V0J,(x) +^Ee^*^(TV^ + ^sin^V(/.)J±(x) 



(5) 



where Jq;(x) {a — ±,z) are the spin currents of the electron; 

Ja(x,r) = - — [(aV«Va) - (VaV^a)], (6) 

with (T± = (Tj; ± icr^;. 

The total Lagrangean written in the local frame is therefore given by 

^ = E «L(5r + ek.)ak, + ^ ^ M2 + (7) 
where the electron energy is now spin dependent, ek± = (k^/2m) =F t^-^ — ^f- 



B. Effective action 

The effective action of the magnetization is obtained by integrating out the electron 
degree of freedom treating the interaction term Hi^t perturbatively. This is reasonable since 
the variation of the magnetization vector is generally very slow compared to that of electron. 
For the case of a domain wall with thickness A, this perturbative treatment corresponds to 
the expansion in terms of (kp^X)'^, /cpt being the fermi momentum of the majority spin. 
The effective action for the magnetization is obtained as 

s,s{e, 0, M) = Smf{m) + AS{e, 0, M), (8) 

where Smf is the well-known mean field action of ferromagnet; S'mf = — trln(9^ + eka) + 
PEAU/2)M^, and AS = J dr < Hi^t > -(1/2) J dr J dr' < //i„t(T)//int(T') >. The brack- 
ets indicate the expectation values with respect to the electron. 

In further calculations, the variable M is approximated by its mean field value Mq 
determined by the stationary condition of 5Smf{Mo)/6Mo — 0, i.e., 



The fluctuation SM of M around Mq can be neglected for the study of low energy behavior 
In terms of Mq the fermi momenta of electrons with spin up and down are given as kp± = 
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(2meF)^''^(l ± (f//eF)Mo)^ . The behaviors of quantities Mq and (kp^a) as a function of 
(U/ep) (with the electron number per site, n, being fixed to be unity) are shown in Fig. ^ 



FIG. 4. Mean field solutions of quantities Mq, fepo- = (fepo-a), J = {JM^a/ep), and 
5 = (A;f^ — kp^)/{kp^' + A;f|) as a function of f/ = (U/ep). The ferromagnetism appears for 
f/ > (2/3) and the complete ferromagnetism is realized for U > 1, where kpi vanishes 



as 



The dynamics of {6, (p) is determined by AS, which up to the order of and dr is given 

AS = JdTY: - cose) + ^ (i(V^)2 + (1 - cose)(V</>)^) 

-IJdrJ rfr'^:^ {[(1 - cos^^)V,0] [(1 - cos^')V,0'] < j;(r, x) Ji(r', x') > 

xx' ij 

+ 2e-'^'t'-'^'\Vie -ismeV^(j)){Vje' + I sin eVj(j)') < JV(r,x)ji(r',x') >} , 

(10) 

where n =< (a'''a)x > is the electron number per site and {6',(j)') = (6'(x', r'), (/)(x', r')). 

The correlation functions are calculated by the random phase approximation (RPA) 
to the Coulomb interaction [T^. The Fourier transform of < J+J_ > with the thermal 
frequency uje = 2TTi/P is evaluated, as depicted in Fig. || diagramatically, to be 

2U 



< JV(q,£)Ji(-q, -i) >= C'l- + C\_Ci_j 



' 



(11) 



where 



1 1 

1 1 



Y^{2k + qr{2k + q) 



hi - /k+q, 



fkl - /k+q,T 



Cl_{q,i)^^—Y^{2k + qy 



(12) 



and the irreducible spin susceptibility is given by 



xiMJ)^li: , (13) 

Here /ko- is the fermi distribution function /ko- = l/[e^(''^/^'^~'^F<') + 1], with the fermi energy 
e-pa = ep ± f/Mg. The 2;-component < JzJz > is calculated similarly, but for our purpose 
only its irreducible part is needed, which is given as 



FIG. 5. The diagramatic expression of the spin current correlation function < J_|_J_ > in RPA. 



C. Effective Heisenberg model (local part) 



To make clear the meaning of the effective action, we devide the action into two parts, 
that is local and non-local in time; = A^ioc + A^dis- The local part up to order dr and 

is calculated from the a;^ = 0, g = components of the correlation functions, which are 
obtained as 



lim < JV(q,0)ji(-q,0) > 
lim< j;(q,0)Ji(-q,0)> 



{k 



FT 



607r2m2f/Mn 



n 



m 



6ij. 



(15) 



The local part of the effective action results in [10| 



AS 



loc 



dr / (i X 



cos( 



(V0)^ + sin2^(V0)' 



(16) 



where the spin is 5* = Mo/2 and the exchange coupling or the spin stiffness J ([J/m] in 
MKSA unit) is expressed by the parameters of the itinerant electron as 



J = 



n 



1 - 



307r^mnUMo 



(17) 



The action Eq.(|IB|) is of the same form as that of the ferromagnetic Heisenberg model with 



5* = Mo/2, as has already been indicated |]T3|,|16|. Hence, there is no formal difference from 
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the case of an insulator in the local part of the effective action, but S in the present itinerant 
system is a continuous variable in contrast to the Heisenberg model. The behavior of the 
exchange energy {JM^a/ e-p) is plotted as a function of {U/e-p) in Fig. |^. 



D. Dissipative part (non-local part) 

We now go on to the non-local part AS'dis. This term comes from the second order 
contribution of ifint and it is this term that contains the characteristic feature of the itinerant 
model, the dissipative effect from the conduction electrons. In the following analysis, we 
consider a planar wall with a spin configuration changing only in x-direction in space, and 
the spins lying in the ?/2;-plane; i.e., (p = 7r/2 (see Fig. [l|(a)). In this case, the dissipation 
arises only from < J_|_J_ > term, and then the non-local part is given by 

A^dis = \jdrj dr'^ ^ e^-^(— ')i E l^ni^) " &^{r')\' < J^q, ^) Ji(-q, -i) >, (18) 
where Qqir) = J^x^'^'^^^^- To estimate this effect that dominates at low energy, the ana- 



lytical continuation to the real frequency iuji ^ ± zO is convenient [p^p^ . The important 
contribution is that of the Ohmic dissipation, which is due to the tu-linear term in the imag- 
inary part of the analytically continued correlation function, Im< J(q)J(— q) > \iuji=u;+io as 



u; ^ um. The Ohmic part of AS'dis turns out to be [10 



1 f.^ f.^> 1 V- I0q(^) - 0q(^')P /Im < JUq)JU-q) > 



as 



A5dis = 7 /rfr/rf/4V '^q\^ 7: ^' lim r^^^"" + ^^^"-^ h;^ic.+.u . ^^g^ 
4 7 7 V ^ (r - r')2 '^^o \^ ncu J 

After straightfoward algebra, the Ohmic contribution to the imaginary part is calculated 

m 



Im < J^(q)J^-q) > |^+,o = re J 7~ ,^3^^ gst(g) + 0{uj'). (20) 
where the function 6st{q), 

^ , , ,1 (^FT - kpi) < \q\ < (A;ft + kpi) 

Ost{q) = S , (21) 

otherwise 
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indicates that the Ohmic dissipation is due to the Stoner excitation, which is a gapless 
excitation that flips a spin of an electron at the fermi surface. The effect of dissipation at 
zero temperature is finally given as 

^"^dis = ' „2 / T^TZi^^M dr dr , J, . (22) 



327r2 J (27r)3 |g|3 ^^^^V J (r - r')^ 

The total action for the direction of the magnetisation vector is AS\oc + A^dis- Since ASdis 
is positive definite, the tunneling rate is seen to be always reduced by this dissipation effect 
within the present semi-classical approximation. In the case of a complete ferromagnet, 
(U/ep) > 1, kpi = and then A^dis is vanishing. 

III. EFFECTS OF DISSIPATION ON MQT OF A DOMAIN WALL 

In the previous section, we have derived based on the Hubbard model an effective action 
A^ioc + ASdis describing the slowly varying direction of the magnetization vector. By use of 
this effective action, we will investigate the quantum tunneling of a domain wall in the case 
as shown in Fig.|I|(a). To do this, we need first to describe the domain wall motion in terms 
of the tunneling variable, the coordinate of the the wall Q. The dissipative effect A^dis is 
to be estimated by use of the variable Q. 



A. MQT of a domain wall 



Let us first study the local part of the action ASioc- Without A^dis, the analysis goes 
the same as in the insulator case 0-^. For the domain wall to exist, anisotropy energy is 
needed, which we shall add phenomenologically to the action AS'ioc as the energy — /^[J/m^] 



in the z-direction as an easy axis and +K± in the x-direction as a hard axis [p]9|-[21|l. The 
action we deal with is therefore given by the following (including A^dis) 



s{e, 



dr I d^x 
K 



S ■ 

i—(h(l — cos( 



(V0)2 + sin2 0(V0)^ 



-1152 2 iv^52 i^2^cos2, 
2 2 



+ A5, 



dis- 



(23) 
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The spin, 5* = Mo/2, and the exchange couphng J[J/m] are determined microscopically by 
Eqs-d) and 

Without dissipation, ASdis = 0, the action Eq.(p3D has a classical solution of a planar 
domain wall, where the spin configuration changes only in one direction, which we choose 
depicted in Fig.|I](a). The position of the wall is at x = Q. Because of the 
anisotropy, the spin points to the ±z direction at infinity and it rotates in the yz-plane (due 
to transverse anitsotropy energy that unfavors the x-direction) around x = Q. For a wall 
moving with a small velocity Q <ti c, where c = K±XSa^ /h, the classical solution is given 
by 

X - 0(t) O 

cos ^(x, r) = tanh cos 0(x, r) ^ i— < 1. (24) 

A c 

Here the width of the wall. A, is given by 

A - (25) 

The dynamics of the wall coordinate Q{t) is determined as follows. For the configuration 
of Eg ■ (|2^) , the Lagrangean in Eg . (|23|) (terms in sguare bracket) is written in terms of a 
collective coordinate Q as (adding irrelavant constant) = — (l/2)Mw(5^ where the mass 
of the wall is given by 

M„ . (26) 

Ay^ being the area of the wall. The strength of the pinning potential for Q produced by 
the defect is to be related to the coercive field He of the system. The magnetic field H 
produces a linear potential oc HQ and thus makes the pinning center Q = metastable (see 
Fig. 0(b)). In the experimental situations, the field is set very close to the coercive field He, 
where the energy barrier vanishes, for the tunneling rate to be large enough. In this small 
energy barrier case, the total potential for Q produced by the defect and magnetic field is 
well approximated by a sum of harmonic and cubic potential |^ . 

ViQ) = \m^uIQ' - ^) • (27) 
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Here the attempt frequency uq and the tunnehng barrier width Qo is shown to be expressed 
in terms of macroscopic parameters in the leading order of e = (ifc ^ H) / 1 as H,p3] 



cuo ~ '^ v^e^ ~ 10^^ X yte^ (Hz) ~ 5 x i//^e^ (K) (28) 

Qo = ^v^A, (29) 

where he = Hc/{2h'yS/a^) is the ratio of the coercive field to the magnetic moment per unit 
volume. The typical values of this parameter is — 10~^ ~ 10"^. Note that the barrier 
width Qq is much smaller than the domain wall width A for the case of shallow potential {i.e., 
small e). The barrier height is proportional to the number of spins in the wall = A^X/a^, 
and is expressed as Uu — fio{{h^Sy /a^)Nhce^^ 5 x Nh^S- (K). 

From these considerations, the dynamics of Q is determined by the Lagrangean 

L{Q) = ^M^Q^ + ^M^uj'.Q' (l - ^) • (30) 

The tunneling rate of the variable Q from the metastable state Q = of the potential 
V{Q) is estimated using the bounce solution of L{Q), which describe Q travelling from at 
r = — oo to Qo at r = and goes back to Q = at r = oo P3|. It is explicitly given as 



= (31) 

The exponent of the tunneling rate without dissipation is given by the value of action 
B = J drLiQ) estimated along the bounce solution. Taking into account the prefactor, the 
tunneling rate without dissipation is estimated as |^ 

To = Ae-^, (32) 



where A ~ IQ^^MN^e^ (Hz) and B ~ MNe'^ [g. Because K and e are small, MQT of a 
domain wall is observable even for very large number of spins if the dissipation is negligible. 
For example, for he = 10~^ and = 10^, MQT is observable {e.g., Tq > lO^^Hz) by setting 
the magnetic field within 1% of the coercive field {i.e., e < 10^^) (see dashed line in Fig.^. 
Instead, e < 2 x 10"^ is needed for = 10^. 
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The crossover temperature Tco from the thermal activation to the quantum tunnehng is 
estimated by the relationship B = f/n/^co where Un is the barrier height of the potential. 
For the present problem, the crossover temperature is roughly the same as the attempt 
frequency uq] 

Teo ^ 5 X (K). (33) 



B. Effect of dissipation 

Now we discuss the effect of dissipation due to itinerant electrons, AS'dis- Within the 
calculation based on the bounce solution, the tunneling rate Eq.(|32D is reduced by dissipation 
to be § 

r = Ae-(^+^^^'-) = Toe-'^^'^^ (34) 

Here, to the lowest order in dissipation, the exponent AS'dis is evaluated as the value of the 
action Eg. (^21) estimated along the domain wall solution Eg . (|2^) with Q{t) given by the 
bounce Eg . (^I]) . For the present planar domain wall solution Eg . , 9q is obtained as 

eq(r) = -nA^e-''^^Q^^^—^^^6,,,oS,,,o, (35) 

where g^'s are the i-th component of the momentum q, and AS'dis reads as 



327rAa J J {r — r'Y J(kp^-k-pi)^x x cosh^x' 



where we have used e ^ 1 and the fact that the integral is dominated by small momentum 
transfer; x < 1. In evaluating the time integral, we note that the bounce solution Q{t) is 
localized within |t| < uJq^, and we approximate the function (Q(t) — Qij')) as 



Q{t) - g(r') = { ^ . (37) 

Qo \t\ < ujq \ |r'| > 2uJo ^ 



The action then results in 
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AS-dis = vNe. (38) 



The strength of dissipation is given by 



31n3.,2 ,2 ^2 4^ /•(%T+%i)fA _ 1 1 



ri = (fcpi — kp,) a — dx ^ — , (39) 

IQn ' ^ a J(fcF|-fcF|)f A X cosh a; 

and its asymptotic behaviors are 

7^(A)^:i^(A:i^-fci,)Vx " '^^t-'^n Fi; ^^^^ 



Thus the effect of dissipation is exponentially small for a thick wall A(A;ft ~ ^Fi) ^ 1- In 
the case of an bulk iron, for example, the domain wall thickness is about A ^ 200A and so 
the effect of dissipation from the itinerant electron would be negligible. For a thin wall, on 
the other hand, as realized in e.g., SmCos (A - 12A), 1] may be large. In Fig. y, rj is plotted 
as a function of A/a for ko = {kp^ + kpi)a/2 = 3 and 6 = (fcpt — kpi)a/2kQ = 0.05,0.1 and 
0.2. It is seen that t] can be of the order of 0.1 for A/a ~ 2 and 6 ~ 0.05. In contrast 
to the present metallic case, the Ohmic dissipation vanishes in insulators at absolute zero, 
since in that case, the major source of dissipation is the magnon, which has a gap due to 
the anisotropy. 



FIG. 6. (left) The strength of the dissipation r/ given by Eq.(|39|) as a function of the width of 

the wall (A/a), for 6 = {kp^ - kpi)a/{2ko) = 0.05,0.1 and 0.2 with ko = {kp^ + kpi)a/2 = 3.0. 
FIG. 7. (right) The tunneling rate F (solid line) for a typical case of an metal [t] = 0.1) with 

total number of spins in the wall N = 10^ and 10^, respectively. The dashed lines are the rate Fq 

in the absence of dissipation. The parameter is taken as he = W~^. 

In the presence of dissipation, the tunneling rate F is given by Eq.(P^ with the action 
given by Eg . (^Sf) . For a typical case of metal; = 0.1, F is shown in Fig. |^ for = 10'^ 
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and 10^ with he = 10~^ . The value = 10^ corresponds, for instance, to a material with 
the wall thichness of ~ lOA and the area of 200Ax200A. The rate Tq without dissipation 
defined by Eq.(^), which corresponds to the case of an insulator are also plotted by dashed 
lines. It is seen that to obtain an observable tunneling rate ( > 10~^Hz) in metals, smaller 
value of e by a factor of about 10~^ is needed than in insulators. 



IV. DISSIPATION IN MULTI-BAND MODEL 

Our preceeding calculation shows strong dissipative effects for a thin wall in weak ferro- 
magnet. In convensional bulk metals, these two conditions might not be easy to be satisfied 
simultaneously. However, in realistic situations, the existence of multi-band {e.g., s-d two 
band model) must be taken into account. 

We consider the simplest case of the multi-band, the s-d model, where the localized 
magnetic moment is due to d electron and the current is carried by s electron. Due to the 
s-d mixing, the antiferromagnetic exchange interaction between localized moment M(x) of 
d electrons and the s electron is present PB 



i7,, = ^7 5:M(x)-(cWW(^)).. (41) 

X 

Hence the Lagrangean of the s electron is of the same form as Eq.(^, except that the last 
term in Eq.(|^) is absent here, since the Coulomb interaction among s electron is neglected. 
The calculation is carried out in the same way as in §|| and |T|. Similarly to Eq . (p9D , the 
strength of dissipation due to s elctron is obtained as (^'S'^is = rj^^^Ne) 

^is) = _ (4^))2]v- r;['^';['^^'dx^-^^, (42) 

/ 647r ^ ^ ai(4;'-fcW)§A a; cosh^ x ^ ^ 

where fcpj is the fermi momenta of the s electron. This expression is the same form as rj 
(Eq. (p9|)) p8| , and hence the behavior of r]*^*-* is accordingly read from Fig.|^, and r/'-'^^ can be 



0(0. 1) if the lower bound of the integration, [kp^ — fcpj^)7rA/2, is smaller than unity. Hence if 
the s-d coupling is sufficiently small: {Xkp ^a <^ 1) (a = {gMo/ep '')), large value of r]^'^^ ~ 0.1 
is possible even in strong ferromagnet. The expression of r]^^^ in this limit is 



17 



V 



IGvr V a 



a 



In 



— Xkk^a 
2 ^ 



(43) 



For a = 0.05, (A/a) ~ 2 and /i;p''a = 3, for example, t]*-*-* ~ 0.02. The result indicates that 
in the multi-band systems, dissipation due to Stoner excitation will be significant for a thin 
wall even in the case of strong ferromagnet. 



V. DISSIPATION DUE TO EDDY CURRENT 

Besides the direct coupling to the electron spin current, the moving wall in metals also 
interacts with the charge current (eddy current) via induced electric field governed by Fara- 
day's law as noted by Chudnovsky et al. 0. The dissipation due to this eddy current is 
Ohmic, but this effect turns out to be small for a case of meso- or microscopic wall, as 
discussed below. 

The electric field induced by the change of the magnetization is calculated from the 
Maxwell equation 

V X E = -yUoM, (44) 

where M = (2h'yS/a^)n is the magnetization vector. For the case of a moving domain 
wall, where the angles are given by Eq.(p^, only the y-component Ey is important, which 
is obtained as 

^, . (2n^S\-.^ ftanh^^ \x-Q\<L 
^ ^ ^ [^j^^Sni^-Q) \x-Q\»L 

where L is the linear dimension of the crosssection of the wall {A^ = L"^)- For L ^ X, 
the electric field is almost constant over the distance |x — Ql ^ L and decays for a larger 
distance. 

The electromagnetic coupling to electrons is 

i7EM = e^A-j, (46) 



where the electron current is given by j(x, r) = —{i/2m){c^'S/c — Vc^c), and the vector 
potential is related to the electric field by E = —A. Since E and hence A depend on Q, this 
coupling -f^EM gives rise to dissipation effects on the motion of Q. 

Treating H-^m perturbatively to the second order and integrating out the electron, we 
obtain the following term in the effective action for Q; 

A^ch = -lJdrJ dr'i5]A^(q,r)A4-q,r') <>(q,r)j,(-q,r') >, (47) 

where ^^^(q, t) and j^(q, r) are the Fourier transforms of the vector potential and the current. 
From the behavior of Ey{-x, r), it is easily seen that Ay{q, r) is well approximated as 

^.(q, r) ^ /io f ^) — e^'^-^A^L X 6,^,oK,o (48) 

for Qx ^ L^^ and Ay ^ for ^ L^^. 

We consider the disordered case where L is much larger than the mean free path of 
electron i; 

L:^i, (49) 

since such a case will be of interest in experimental situations. Because of this condition, 
the correlation function < j^(q, T)ji^(— q, r') > can be approximated by its value at g = 0. 
In this case, the Ohmic contribution is written in terms of the conductivity a as 

< jM^ r)j.(-q, r') >^ - H uduje^^^-^'^ = ^ (50) 

TT JO 7r(r — T )^ 

The action is thus written as 



47r 



f2h-fS\ 



aA^L X I dr I rfr'^^^-^^, (51) 



where we have used {Q{t) — Q{t'))/L ^ 1 and subtracted an irrelevant constant. Therefore 
the strength of Ohmic dissipation r]^'^^^ due to the eddy current, defined by AS'ch = r]^'^^^Ne 
{N = Au,X/a^ and e arises from Qo — -\/eA), is ||29[| 
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For a typical metal, a ~ 10^[i7~^m~^], it is •r]^'^'^^ ^ 10"''' x (\/a){L/a) if we choose S* ~ 1 
and a = 3 X 10"^'^ [m]. Hence, as far as (A/a) < 10^, which is commonly satisfied, dissipation 
due to eddy current is small in mesoscopic wall we are interested in; {L/a) < 10^. 



VI. DISCUSSIONS AND CONCLUSION 

We have studied the case of a domain wall where the plane of the wall is parallel to 
the easy (?/z-)plane. It is also possible to consider based on the action Eq.(^) another 
configuration where spin direction changes in the easy {z-) direction of spin as shown in Fig. 
|l|(b) For the MQT of this wall, the results obtained above do not change. 

A contribution from the magnon pole in the correlation function < JJ > leads to 
the dissipation with a gap Aq due to anisotropy. In the present problem of tunnel- 
ing, the effect needs to be evaluated quantitatively p0| , |30[| . The strength of dissipa- 



tion ry(p°^'=), which is defined by the value of the action devided by (A^e), is obtained as 
j^(poie) _ (37r/40)Mo exp(— (Aq/cjo))- Since experiments are usually carried out in highly 
anisotropic materials with Aq/cjo — 10, this contribution is very small compared to the 
Ohmic dissipation for the case of a thin wall. 

Our result shows a distinct difference between MQT of thin walls in metallic and insu- 
lating magnets. Unfortunately the experiments carried out so far appear not yet be able 
to test the theoretical prediction of the effect of dissipation due to itinerant electrons. The 
experimental result on small ferromagnetic particles of Tbo.5Ceo.5Fe2 suggests the motion of 
a domain wall via MQT below Tco — 0.6K [P]. In this experiment, the width of the domain 
wall is about 30A and according to our result, rj oc exp[— 7rA(A;Ft "^f;)], the dissipation from 
electron spin current is negligible for such a thick wall. This may be the reason why the 
result of the crossover temperature ~ 0.6K is roughly in agreement with the theory p| with- 
out dissipation (Eq.(^) with /i^ ~ 4 x 10~^ and e ~ 0.1). On the other hand, the domain 
wall in SmCo5 is very thin A ^ 12A, and our result suggests strong effect of dissipation due 
to Stoner excitation, which will be interesting to observe. Experiments on bulk crystal of 
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SmCos.sCui.s with very thin walls (a few times a) have been performed 0, although quan- 
titative comparison is not easy since many walls will participate in the relaxation processes 
of the magnetization in these experiments. 

To conclude, we have studied the macroscopic quantum tunneling of a domain wall in a 
metallic ferromagnet based on the Hubbard model. The effective action of the magnetization 
is obtained to be the same form as a ferromagnetic Heisenberg model but with a non-local 
(retarded) part, that describes the dissipative effect on magnetization. In contrast to the 
case of insulators, the Ohmic dissipation exists even at zero temperature because of the 
gapless Stoner excitation in the itinerant electron. The effect is negligible for a thick domain 
wall where experiments so far have been carried out. On the other hand, important effects 
of the Ohmic dissipation are expected in thin domain walls with thickness A comparable to 
the inverse of the difference of the fermi momenta (fcpi — k-pi)~^. In the strong ferromagnets, 
the Ohmic dissipation arising from the magnetic band is weakend. In this case the existence 
of multi-band needs to be taken into account whose effect can be significant even in strong 
ferromagnets. Dissipation due to the eddy current has been shown to be negligibly small 
for a mesoscopic system. We believe the observation of the dissipation due to the Stoner 
excitation is within the present experimental attainability. 

In the case of a single domain (ferro-)magnets, where another typical MQT can occur 



| 31| , our analysis implies that the dissipation due to the conduction electron is irrelevant 
since only the weak component of the electron excitation {i.e., q = 0) is coupled to the 
uniform magnetization. This fact is seen as a vanishing rj (Eq. (0)) in the limit of single 
domain magnet, A ^ oo. 

We have investigated the MQT in metal paying attention to the role of the conduction 
electron. On the other hand the electron transport properties should be affected by the 
MQP of magnetism. This new possibility will make the qunatum tunneling in magnetic 
system more exciting in metals. In fact the experimental study in this direction has been 
started recently [^|, where the resistivity of a mesoscopic Ni wire was measured in the 



magnetic field. Small jumps have been observed there, which may be related to the MQT 
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of domain walls at low temperature. The novel feature of this experiment is that it is on a 
single mesoscopic sample, hence one will be able to detect a single MQP event. This is in 
contrast to the magnetization measurement using SQUID where in general a large number 
of small particles are needed to get a sufficient signal, and hence one must worry about the 
distribution of energy barriers in its interpretation. Another new possiblity in the study of 
the mutual effects between the conduction electron and the MQP of magnetism would be 



magnetic dots in semiconductors |]33|, where one can control the interaction by changing the 
electron density. 

The MQP in magnetic metals and semiconductors will be a very important field also 
from the point of view of application and there are, we believe, a lot of work to be done in 
this new field. 
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